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BIOELECTRIC SIGNAL ANALYSIS AND MEASUREMENT 

The goals of this research project are: (a) to use nonstationary 

time-series techniques to analyze EEG signals for the estimation of 
alertness; (b) to extract time-varying order in sequential time-series 
measurement of these data; and (c) to devise strategies for obtaining 
optimal representation of the EEG signal. Significant accomplishments 
on each of these, goals have been made. Most of the results have been 
presented an various international and national scientific conferences 
a-a published in the open literature. This report will mainly consist 
of reprints of these published papers. However, we shall give a brief 
account of our accomplishments as follows. 

<i) Vie have successfully implemented a frequency-discrimination technique 
tor EEC- to obtain the rate of change of its phase. This is based' on the 
m-ccei that EEG is a phase and amplitude modulated signal with carrier 
frequency center around the or-rhythm, viz., 

e(t) = a(t) cos ^2TTf^t + <f)(t)J 

where a(t) represents the slow amplitude variation of the signal; 4>(t) 
represents the phase variation of the signal; and f designates the 
a -rhythm frequency which is about 8 - 13 Hz, Details are reported in 
a paper published in the Proceedings of the 1972 National Electronics 
Conference. 

(2) A mathematical model for the photically stimulated EEG signals has 
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been conceived and simulated, based on the phase alignment of EEG signals 

with the photic stimulus. In the model, a parameter b which represents 

\ 

the speed of the propagation of the phase uncertainty is thought to be 
related to the alertness states * Calibration of this parameter with 
respect to the alertness stages would make this parameter a useful tool 
for alertness estimate. Details of this model are described in a paper 
entitled -•* Model zor the Photicallv Stimulated Electroencephalographic 
Signals'', published in the Proceedings of the 12 th Annual San Diego 
Biomedical Symposium. 


(o) Spectral analysis of LilG signals has been used for many years* 
However, real rime nard copy displays of the time-varying power spectrum 
that permit tracking of dynamic brain states have been lacking. A 
scheme was implemented. This scheme enables us to obtain a display with 


the amplitude being exhibited both in waveform and in' gray intensity. 
Visually, this hard copy produces a three-dimensional effect. The display 
shows clearly the dynamic changes in the power spectrum of the EEG* 

The detailed description of this program dubbed n GIFBUF n written by 
F. Mansfield . was submitted earlier. 


(4) The representation of EEG signals has always been a major problem 
in the automatic analysis of EEG signals. We have developed a repre- 
sentation for EEG signals and their in-phase and quadrature components* 

This representation makes easy the monitoring and tracking of EEG frequencies 
and phases. The development and implementation of this representation 
resulted in two papers : "Real-time EEG Analysis and Monitoring Using’ 
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In-phase and Quadrature Components", published in the Proceedings of the 
26 1 Annual Conference on Engineering in Medicine and Biology, and 
"Error-free Representation of EEG Signals", published in the Proceedings 
of the 1973 IEEE International Conference on Systems, Man and Cybernetics. 


Co) It has been a major task to separate signal from noise in EEG analysis 
** devised a scheme by modeling an EEG signal in the form of a sinusoid 
waveform and an additive noise which represents the incoherent component 
a ~ vne 3i S r --^- Wrth this model, we developed a scheme which will estimate 
•_ne signal— to— noise ratio in the most commonly used processing situation; 
Tiz., averaging. The technique and results are described in "Estimating 
Signal and Noise in Coherent Time Averages' of EEG Data", published in 
me Proceedings of the 26 Annual Conference oh Engineering in Medicine 


ioj A nonlinear oscillator model for EEG signals has been developed. 

This model has proven to encompass many reported phenomena and predicted 
several unreported phenomena when the subjects were under periodic 
photic stimulation. The oscillator used is the van der Pol oscillator 
described by: 

/ 2 \ 90 

x - \i\l - x Jx + w Q x = oj Q e(t) 

where x(t) denotes the EEG signal to be modeled; is the autonomous ' 

alpha frequency; e(t) is the external stimulus; and p is the coupling 
coefficient. This model also elucidates the effect of a stimulus flash 
to the phase of an EEG signal. The details are reported in a paper, 

"A Nonlinear Model of EEG Entrainment by Periodic Photic Stimulation", 
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published in the Proceedings of the 7 th Annual Conference of the 
Neuroelectric Society , 

(7) To monitor or track the alertness states via EEG signals, we need 
the ability to predict the EEG waveforsm. For this purpose, we used 
an autoregressive process representation. This endeavor resulted in a 
paper entitled, "Prediction of EEG Waveforms by Using an Autoregressive 
Model , lo be presented at the 1975 San Deigo Biomedical Symposium* 

Tne lull paper will be published in the Proceedings * A reprint of the 
abstract of this paper is attached. 

In summary, we have accomplished most ox the objectives as proposed* 


4 



APPLICATION OF FREQUENCY DISCRIMINATION TECHNIQUE TO THE 
ANALYSIS OF ELECTROENCEPHALOGRAPHIC SIGNALS 


D, C. Lai 

Stanford University 
Department of Electrical Engineering 
Stanford, California 94305 

R. L. Lux 

University of Utah Medical Center 
Cardiology Division 
Salt Lake City, Utah 84112 


Reprinted from the Proceedings of the National 
Electronics Conference, vol. 27, pp. 80-85, 
October 1972. 


A portion of this work was carried out at the Department 
of Electrical Engineering, the University of Vermont, Bur- 
lington, Vermont, and was supported in part by NASA Grant 
NGR 46-001-038, At Stanford University, it was also sup- 
ported in part by NASA Grant NGR 05-020-575 from NASA -Ames 
Research Center and ARPA Contract DAHC15-72-C-0232. 



APPLICATION OF FREQUENCY DISCRIMINATION TECHNIQUE TO THE 


ANALYSIS OF ELECTROENCEPHALOGRAPH! C SIGNALS* 


D. C. Lai 


R . L . Lux 


Stanford University 

Department of Electrical Engineering 
Stanford , California 94305 


University of Utah Medical Center 

Cardiology Division 
Salt Lake City, Utah 84112 


INDEX TERMS - EEG , Frequency discrimination , Narrow - 
band process , Entrainment of EEG alpha rhythm , Com - 
puter analysis. 


I . INTRODUCTION 

The experimental analysis of EEG waveforms has 
been dominated by spectral decomposition and ampli- 
tude analysis, mainly autocorreiogram techniques 
though crosscorrelation technique has been used for 
phase measurements. These techniques as applied to 
the analysis of EEG signals are abundant in litera- 
ture. To cite a few of thorn, we note that the ear- 
liest ones were done by Brazier and Barlow [l], [2], 
[3], [4], and notable simplifications done by Kamp, 
et al . [5], DeBoer and Kuyper [6], and Vo-Ngo^ : et 
al [7], The fact that the phases of EEG signals as 
compared to a reference clock signal contain infor- 
mation about the mental alertness has been pointed 
out to us by Anlikcr [8] who has used phase-vector 
and contour-graphic techniques to study the phase 
entrainment phenomenon. Adey and Walter [9] have 
used a phase detection technique in the analysis of 
EEG records in the cat. The EEG a-rhythm has been 
treated as the result of mutual synchronization of 
a population of spontaneously oscillatory processes 
by Wiener [10] , The potential usefulness of the 
phase entrainment in either recognizing the patho- 
logical states of the brain or classifying mental 
states might be widely explored if investigators had 
a reliably simple direct measurement technique. In 
this paper' we describe a frequency discrimination 
technique as realized by an inexpensive frequency 
discriminator f or direct on-line measurement of the 
entrainment phenomenon in EEG as entrained by the 
frequency of a sensory stimulus. We demonstrate the 
use of this device for detecting the presence or ab- 


carrier f requency . The process <Kt) can also 

be interpreted as a relative phase angle, in the 
sense that the alpha signal process a(t) differs 
in phase from the signal (cos 2-n^t) by <t>(t) . It 
is only in this manner that we are able to ascribe 
meaning to the phase of EEG as a time-varying quan- 
tity. The fact that the angular velocity or fre- 
quency in radians per second is the time derivative 
of the angular position leads us to define the in- 
stantaneous frequency f^(t) in cycles per second 
or Hertz by 


f. (t) 
x 


_L it 

2tt dt 


( 2 ) 


The second term in Equation (2) which is proportion- 
al to the rate of change of the time-varying phase 
4>(t) can be interpreted as the instantaneous fre- 
quency deviation relative to the a frequency 
provided it is stable. Our technique to be describ- 
ed in detail provides a direct measurement of (d$)/ 
(dt) . The study on the statistical properties of 
the random processes a (t) and Mt) will be re- 
ported in another article. 

When the brain is excited by an outside stimu- 
lus frequency f Q and the EEG is entrained by the 
stimulus, then we should have 


f . (t) 


f a = 


t + T < t < t r 
n n — — f 


(3) 


where t n and tf signify respectively the time 
instants at which the stimulus is on and the stimu- 
lus is off, and r n denotes the time delay between 
t n and the time instant at which the entrainment 
occurs and Tf, the time delay between tf and the 
time instant at which the entrainment disappears . 
Equation (3) implies that 


sencc of the stimulus effect and the measurement of 
the time delays in the entrainment . 

In our study, the EEG signal is first filtered 
by a narrow-band filter with center frequency about 
the Ci-rhythm of the individual , We shall denote 
this filtered EEG signal as alpha signal a(t). It 
is then reasonable to consider the alpha signal as 
narrow-band random process with the alpha frequency 
% as the mean frequency of the spectral band. Then 
a sample function of this random process is express- 
ed as [11] 

Ct(t) = a (t) cos [2TTf a t+0(t>] (1) 

If the bandwidth of its power density spectrum is 
much smaller than its mean frequency then the 

processes a(t) and Mt) in Equation (1) will be 
slowly varying functions of time as compared to 
cos 2TTf a t so that the interpretation of a(t) and 
<t>(t) as envelope and phase has meaning. The alpha 
signal process rt(t) can thus be interpreted as 
amplitude and phase modulated signal process with 


d<t> 

dt 


= o, 


+ T n < 4 S *1 


+ T, 


(4) 


Since our technique measures (d<l>)/(dt) directly, 
the measurement of EEG signals which were recorded 
while the human subjects with eyes closed were stim- 
ulated by stroboscopic flashes at the rate of 10 
flashes per second for one minute, then no stimula- 
tion for one minute, then another minute of stimula- 
tion, etc., does show that (d<J>)/(dt) = 0 for al- 
ternate one minute intervals. 

The schemes used for the measurement of (d4>)/ 
(dt) , T n , and Tf , and for the automatic detection 
of the states of stimulus-on and stimulus-off in the 
EEG are to be described next. 


II . METHOD 

An ideal frequency discriminator should pro- 
duce an output voltage linearly dependent on input 


$ 
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frequency. There are many ways for the realization 
of frequency <1 is crimination. The most commonly used 
one is the so-called balanced demodulator with the 
well-known S curve. The scheme which we have used 
here is the zero-crossing detector type as described 
in [12]. In order to eliminate the influence of the 
amplitude variation a(t), we use an amplifier fol- 
lowed by a hard limiter in front of the detector. 

The output of the frequency discriminator is propor- 
tional to -the instantaneous frequency fi<t) as ex- 
pressed in Equation (2), and does not depend on the 
amplitude of the signal. We are, however, interest- 
ed in the rate of change of 4>(t). The value % 
in Equation (2) is subtracted. In our device ifo 
can be preset at any value in the range of 5-15 Hz. 
The device thus yields (K/2ir) (d$/dt) on-line. The 
characteristic curve of the device is given in Fig.l. 
In Fig. 2, we show the functional block diagram of 
the analysis scheme. Detailed circuitry of the de- 
vice is given in Fig. 3. 

The response time or delay T n (or if) is 
defined as the length of time elapsed since the on- 
set (or end) of stimulus to the time instant at 
which the brain wave is entrained (or desynchroniz- 
ed) , This delay T n (or Tf) is, in a sense, the 
time 1 which the brain takes to synchronize (or desyn- 
chronize) . There is a noticeable delay between the 
onset of stimulus-on (or stimulus-off) and the time 
at which the entrainment is measurable. This fact 
is shown in Fig. 4 where the (d4>)/(dt) as the out- 
put of the frequency discriminator does not reach a 
flat plateau (or start to wander) immediately after 
the onset of stimulus-on (or stimulus-off) . The 
fluctuation of the x(t) £ (db)/(dt) curve is mea- 
sured by the quantity 

P = ^ |x(t > - m l (5) 

1=1 


within a specified window width N, where m N 
the mean value of (db)/(dt) or x(t) in the 
window; i.e., 


m 


N 



x(t i ) 


is 


( 6 ) 


New values of p are calculated by sliding the win- 
dow. The onset of the flatness (or the wandering) 
is determined by the time instant at which p is 
below (or above) a preset threshold that is deter- 
mined experimentally. The response time or delay 
T n (or x f ) is measured as the time lapse between 
the onset of stimulus-on (or stimulus-off) and the 
onset of the flatness (or the wandering) of (d<»)/ 
(dt) as measured by p. Owing to the delay intro- 
duced by the window width, a correction constant is 
subtracted from the above obtained time intervals. 

The detection scheme utilizes the same program 
as that used in measuring the response time or de- 
lay, but without the prior knowledge of the onset 
of stimulus-on and stimulus-off. In other words, 
the influence of stimulus (entrainment) is said to 
be detected whenever the quantity p reaches a val- 
ue below a preset threshold. 

These schemes are realized by digital computer 
programs in PAL III language used on a PDP-8 compu- 
ter. Both of these measurements can be performed 
on-line. Figure 5 depicts the flow chart of the 
computer program for the response— time measurement. 

III. RESULTS 

Although the measurements were performed on EEG 
signals recorded on magnetic tape, they can be per- 


formed on-line. As mentioned before, these EEG sig- 
nals were recorded while the subjects with eyes 
closed were stimulated by stroboscopic flashes for 
one minute at the rate of 10 flashes per second, 
then no stimulation for one minute, then another min- 
ute of stimulation, etc. The flashes were generated 
by Grass photo stimulator model PS-2E on its inten- 
sity scale No, 2. The EEG signals are the differen- 
tial potentials between the electrodes placed at 
left parietal and left occipital. The EEG signals 
are referred to as either stimulus-entrained or non- 
entrained. In order that the EEG signal be more ap- 
propriately described by the narrow-band-process mod- 
el, the EEG signal was passed through a narrow-band 
filter with a band-width of 1 . 5 Hz centered at 10 Hz 
to obtain the alpha signal a(t) as described pre- 
viously. 

A section of the output of the frequency dis- 
criminator along with the same section of EEG sig- 
nals of subject B being analyzed and a square wave 
indicating when the stimulus was turned on or off 
are shown in Figs. 6 and 7. The display of apparent 
high frequency components in the raw EEG of Figs. 6 
and 7 are due to the slower writing speed (30 cm per 
minute) of our strip chart recorder as compared to 
the conventional writing speed. It is clearly seen 
that there is a delay following the onset of the 
stimulus before entrainment is detected; the EEG re- 
mains synchronized for a short period of time also 
after the stimulus was turned off. In Figs. 8, 9, 
and 10, we show the output of the frequency discrim- 
inator with the early, middle, and late sections of 
EEG signals of subject H as the input along with a 
square wave indicating the events of stimulus-on and 
stimulus-off. The offset of time coordinates in the 
above figures is due to the positions of pens of the 
strip chart recorder. From these figures, we observe 
that subject H synchronized with the stimulus very 
well at the beginning section (13th min. to 19th 
min.) of the experiment as shown by the measurement 
of (d4>)/(dt) depicted in Fig. 8. As shown in Fig. 
9, however, the subject T s EEG signal did not entrain 
to the stimulus frequency as well during the middle 
(105th min. to 111th min.) of an experimental ses- 
sion which lasted about three hours. There is no 
sign of entrainment of the EEG signal in the end 
section (176th min. to 182nd min.) of the experiment 
as shown in Fig. 10. The capability to synchronize 
with the stimulus may be a measure of certain brain 
states such as alertness. Further study is needed 
for this measurement by this technique. The fre- 
quency discrimination technique as realized by an 
inexpensive device is effective in measuring the en- 
trainment phenomenon of EEG signals. 

It is possible that variations in the delays 
T n and Tf in synchronization and desynchroniza- 
tion, respectively, might give some indication of 
change of levels of alertness, although our data 
does not include alertness estimates. Results of 
the measurement of delays for both subjects are 
shown in Figs. 11 and 12, respectively. The abscis- 
sas in the figures denote the time since the begin- 
ning of the experiments and the ordinates are the 
measured delays. Any x n > 25 seconds is considered 
as not synchronized. Subject B seems to have longer 
response time to stimulus-on in the first hour of 
experiment than those in the second hour. Subject H 
did not entrain so well after the first hour of ex- 
periment. In general, we see that it takes longer 
for the subjects to synchronize with the stimuli 
than it takes to return to the natural states after 
the stimuli were shut off. The statistics of x n 
and Tf and their correlation with the mental 
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states of subjects need farther study. However, the 
technique for measuring 'these delays is satisfactory. 

In Figure 13, we show a typical section of the 
conditions or states of the stimulus-on and stimu- 
lus-off of EEG as results of computer analysis of 
the output of the frequency discriminator in com- 
parison with the timing track for the on and off of 
the stimulus. This was generated by a double- 
threshold detector. The thresholds were determined 
experimentally. We analyzed the records of both 
subject B and subject H. The ratio of the sum of 
the differences between the time of the on and off 
conditions of EEG as determined by the detector and 
on-off timing track for the stimulus to the total 
time of on and off as measured by the timing track 
was calculated. This ratio is 28,6 per cent for 
subject B and 35.8 per cent for subject H. We did 
not take into account the time delays introduced by 
the smoothing window and the frequency discrimina- 
tor in the evaluation of these ratios. The values 
of these ratios may be reduced. The use of the 
above-mentioned ratio as a performance measure of 
the detector is appropriate if we do not have the 
knowledge of the length of time and the manner in 
which the stimulus is turned on or off. With the 
knowledge that the stimulus was turned on for one 
minute, then shut off for one minute, then was on 
for another minute, etc., and counting only the 
number of the on-off cycles as determined by the 
detector, we obtain the following result: 

(1) For subject B, there are 57 on-off cycles detec- 
ted by the detector as compared to 57 on-off cycles 
on the timing track (a 100 per cent detection) , 

(2) For subject H r there are 56 on-off cycles detec- 
ted by the detector as compared to 59 on-off cycles 
on the timing track (a 94,7 per cent detection). 

IV. CONCLUSIONS 

The frequency discrimination technique describ- 
ed here has been effective in measuring the entrain- 
ment phenomenon in EEG signals. We have shown that 
the narrow-band-process model used here for charac- 
terizing signals is valid. The technique as real- 
ized by an inexpensive device described previously 
provides an effective way for further study on the 
synchronization phenomenon as related to alertness. 
Statistical properties of the phase <Kt) must be 
gathered for further interpretation of mental states. 


V. ACKNOWLEDGEMENT 

The authors wish to thank Dr. J. E. Anliker of 
NASA- Ames Research Center for furnishing the EEG 
data used in this study, and for many valuable 
suggestions . 


REFERENCES 

1. M.A.B. Brazier, Crosscorrelation and Autocorre- 
lation Studies of Electroencephalographic Poten 
tials, Electroenceph . Clin. Neurophysiol . , 1952 
4:201. 

2. M.A.B. Brazier and J.S. Barlow, Correlation 
Analysis of Brain Potentials, Electroenceph. 
Clin. Neurophysiol., 1955, 7:463. 

3. J.S. Barlow, et al . , The Application of Auto- 
correlation Analysis to Electroencephalography, 
Proc . of the First National Biophysics Conf , , 
Yale Univ, Press, New Haven, 1959, 622p. 

4. J.S. Barlow, Auto-correlation and Cross-corre- 
lation Analysis in Electroencephalography, IRE 
Trans, on Med. Electron,, 1959, ME-6; 179-183. 

5. A. Kamp, W. Storm Van Leeuwen and A.M. Tielen, 

A Method for Auto- and Cross-relation Analysis 
of the EEG, Electroenceph. Clin. Neurophysiol., 
1965, 19:91-95. 

6. E. DeBaer and P. Kuyper, Triggered Correlation, 

, IEEE Trans, on Bio-Med . Engineering, 1968, BME- 
1 15:169-179. 

7. B. Vo-Ngoc, et al., A Possible Procedure to Cal 
culate Correlation Functions for the EEG (Cor- 
respondence) , IEEE Trans, on Bio-Med. Engineer- 
ing, 1970, BME-17: 265-267. 

8. J.E. Anliker, Private Communication 1968. 

9. W.R. Adey and D.O. Walter, Application of Phase 
Detection and Averaging Techniques in Computer 
Analysis of EEG Records in the Cat, Experimen- 
tal Neurology, 1963, 7:186-209. 

10. N. Wiener, Nonlinear Problems in Random Theory , 
M.I.T. Press, Cambridge, Mass., 1958, 67-77. 

11. W.B. Davenport , Jr . , and W.L. Rooty An Intro . 
to the Theory of Random Signals and Noise , 
McGraw Hill, N.Y. , 1958, 158 p. 

12. M. Schwartz, Information Transmission, Modula- 
tion, and Noise , Second Ed., McGraw Hill, N.Y. , 
1970, 261-263. 


82 




VOLTS 

FIG. I FREQUENCY DISCRIMINATOR OUTPUT CHARACTERISTIC. 



FIG. 3 FREQUENCY DISCRIMINATOR - CIRCUIT DIAGRAM. 



FIG. 4 RESPONSE TIME OR DELAY MEASURED. 



FIG. 5 FLOW CHART FOR 
THE RESPONSE-TIME MEASUREMENT PROGRAM. 
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PIG 6 MEASUREMENT OF — AND THE DISPLAY OF RAW FIG. 9 MEASUREMENT OF ~ OF THE ALPHA SIGNAL OF THE 

dt ax 

EEG INPUT TO THE SCHEMATIC DIAGRAM SHOWN IN FIG. 2 MIDDLE PART OF EEG RECORD OF SUBJECT H. 

WHEN STIMULUS IS ON. STRIP-CHART RECORDER WRITING 
SPEED IS 30 CM PER MINUTE. 



FIG. 7 MEASUREMENT 1 OF — AND THE DISPLAY OF RAW 
dt 

EEG INPUT TO THE SCHEMATIC DIAGRAM SHOWN IN FIG. 2 
WHEN STIMULUS IS OFF. STRIP-CHART RECORDER WRITING 
SPEED IS 30 CM PER MINUTE. 


FIG. 10 MEASUREMENT- OF ^ OF- THE ALPHA SIGNAL. OF 
dt 

THE LATE PART OF EEG RECORD OF SUBJECT H. 




FIG. 8 MEASUREMENT OF -± OF THE ALPHA SIGNAL OF THE 
dt 

EARLY PART OF EEG RECORD OF SUBJECT H . 


0 



TIMING TRACK 



FIG. 13 THE STATES OF STIMULUS -ON AND STIMULUS -OFF 
OF EEG AS RESULTS OF THE ANALYSIS IN COMPARISON 
WITH THE ACTUAL ON AND OFF OF THE STIMULUS . 
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Fig. 11. RESPONSE TIME STUDY OF SUBJECT B. 



Fig. 12. RESPONSE TIME STUDY OF SUBJECT H. 
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REAL-TIME EEG ANALYSIS AND MONITORING 
USING IN-PHASE AND QUADRATURE COMPONENTS 


44.5 


M. Ein-Gal and D.C. Lai 


In the study of brain functions or monitoring 
patients in clinical laboratory, a non-invasive 
scheme for monitoring either psycho-physiological 
states or pathological states via the analysis of 
EEG signals is an important tool. In view of the 
time-varying nature of these states, we are inter- 
ested in monitoring them on a moment -by-moment ba- 
sis. To this aim, we have devised a digital scheme 
which is capable of tracking the amplitude, fre- 
quency, and phase of a particular EEG rhythm in 
real time. The amplitude, frequency, and phase are 
assumed to be low-pass process. These 3-tuples may 
be treated as a 3-dimensional vector whose trajec- 
tories in the 3-space relate to the change of 
psycho-physiological states or pathological states. 
This scheme is implemented on a PDP- 15 computer system* 
We characterize EEG e(t)as a narrowband process, hence 

e (t ) = a (t ) cos 0(t) = a(t) cos [u Q t + cp (t ) ] , 

or e(t) = c(t) cos go t + s(t) sin go t 
o o 

By tracking the center frequency go^ , we may de- 
compose the signal into two components; viz. , the 
in-phase component c(t) and the quadrature com- 
ponent s(t) . A digital system depicted in Fig. 

1 tracks the center frequency go and resolves the 
input EEG waveform into c(t) and s(t) . This 
system is composed of three major components: 

(a) Voltage-controlled oscillators and 
amplitude regulator, 

(b) Frequency discriminator, 
and (c) Low pass filters. 

The VCO is realized by a second-order difference 
equation. The locations of the poles are controlled 
by the amplitude regulator and by the frequency 
discriminator. The realization of the frequency 
discriminator is obtained by averaging the rate of 
change of phase over a preselected interval in 
accordance with the desired time constant of the 
loop. Applying c(t) and s(t) to the horizontal 
and vertical axes respectively, we obtain a two- 
dimensional trajectorial plot where the envelope 
a(t) and the phase Cp(t) are proportional to the 
radius and angle respectively in the polar coordin- 
ate system. 

As an example, this real-time monitoring tech- 
nique was applied to the EEG data obtained while 
the subjects with eyes closed were stimulated by 
stroboscopic flashes for one minute at the rate of 
10 flashes per second, then no stimulation for one 
minute, then another minute of stimulation, etc. 

The flashes were generated by Grass photo stimula- 
tor model PS-2E on its intensity scale no. 2. The 
EEG signals are the differential potentials be- 
tween the electrodes placed at left parietal and 
left occipital. The results are shown in Figures 2 
and 3. Figure 2 depicts the representation of EEG 
when stimulus was off by using this scheme. Figure 
3 shows the representation of EEG when stimulus was 
on. The plots in both figures are normalized with 
respect to certain scale for visual convenience and 
they are plotted simultaneously. The top one is 
the trajectorial plot of the envelope a(t) and 
phase Cp(t) in polar coordinates. The second one 
is the plot of frequency versus time. The third 


and forth are the rate of change of amplitude and 
phase, respectively. The fifth and the sixth are 
the phase Cp(t) and amplitude a(t) , respectively, 
The seventh is the EEG signal analyzed. The last 
is the output of VCO which is used as reference. 
Notice that the frequency and the phase are practi- 
cally constant when stimulus was on. 



Figure 1. Resolver of 
in-phase and quadrature signals 



Figure 2. Representation of 
stimulus-off EEG signal 



Figure 3. Representation of 
stimulus-on EEG signal 
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ERROR-FREE EEG SIGNAL REPRESENTATION 
M. Ein-Gal and D. C. Lai 
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Summary 

The representation of elect roencephalographic 
signals is an important problem for real-time moni- 
toring of either psycho-physiological states or patho- 
logical states. A representation of any EEG rhythm 
in its in-phase and quadrature components is given. 

This representation gives zero error when the signal 
is reconstituted from its representatives. Deriva- 
tion of the scheme is reported. Using this scheme, a 
digital real-time system for obtaining this represen- 
tation is designed and realized by a PDP-15 computer. 
An example with a typical EEG signal as the input 
to the system is included to demonstrate the accuracy 
of the representation. 

I. Introduction 

A digital real-time system is described for 
obtaining the representation of a particular EEG 
rhythm in its in-phase and quadrature components with 
the constraint that the error is zero when synthesized. 
In a way, this system provides a form of data com- 
pression. The resultant representation can be advan- 
tageously employed for data transmission and data 
storage as well as real-time monitoring Of either 
psycho-physiological states or pathological states. 

For these applications, the requirement of an error- 
free synthesis or reproduction of the EEG signal is 
imperative. 


We may write 

y(N) = X T (N)H(N) (2) 

where X(N) is a two-dimensional vector whose compo- 
nents x^ (N ) and x^CN) are the slowly-varying compo- 
nents sought after. Note that the Equation (2) is not 
unique since there exist many X(N) which could satisfy 
(2). For certain particular choice, x^N) and x (N) 
are related by Hilbert transform. Let X(N|N-1) Be a 
function of the past values of X(N). For instance, 
X(N|N-1) could be a linear combination of X(0), X(l), 

. X(N-1). Denote the difference between X(N) and 
its X(N|N-1) by 

E (N) = X (N) - X<n|N-1) . (3) 

The norm of the difference is 

||e(N)|| = [VfNmioJ 14 . (4) 

We choose X(N) by minimizing (4) but satisfying (2) 
simultaneously. The solution thus obtained is 

X 00 = 1 - A(N)j X(n|n- 1) + H(N)y(N) (5) 

where I is the identity matrix and A(n) is the matrix 
with its components given by the double-f requency . 
oscillator; i. e. , 


A common scheme for obtaining the in-phase and 
quadrature components of a signal is to multiply 
the signal by an oscillator outputs at quadrature with 
the center frequency f Q of the signal and then passing 
the resultants through low-pass filters in the for- 
ward path. However, this method does not provide a 
faithful replica of the signal when the two compo- 
nents are synthesized. The error arises mainly from 
the phase distortions of the filter or the pure delays 
of the non-recursive filter. Our digital system uses, 
instead, an oscillator running at twice the center 
frequency of the signal and a filter in the feedback 
path. This filter does not affect the reproduction 
fidelity; however, it plays a role in determining 
certain statistical properties of the slowly varying 
components of the resultant representatives. Hence, 
one may choose a filter to minimize certain proper- 
ties of the outcome such as the bandwidth, variance, 
etc. 


II. Derivation of the Scheme 

Let the input signal in its digitized form be 
y(N), We shall define a two-dimensional oscillator 
vector H(N) by 

H T (N) = [coe2Tcf 


This work w'as supported in part by the Advanced 
Research Projects Agency of the Department of Defense 
under Contract DAHC 15-72-C-0232 and NASA under 
Grant NGR— 05— 020— 575. 
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If X(n|n- 1) is the best linear estimator of X(N) based 
on its past history, then the representation of the 
input signal y<N) by the vector E(N) corresponds to 
the innovation representation v(N) with 


||E(N)|| 2 = v 2 (N) (7) 

where y(N) = y(N) - $(n|n- 1) and y(N|N-l) is the best 
linear estimator of y(N) in terms of its past values 
y (0) , y (l), . . . , y(N-l). 


III. Results and Conclusions 


The digital system is depicted by a block diagram 
shown in Figure 1. The oscillator is realized by a 
second order difference equation whose coefficients 
are controlled by the amplitude regulator and the 
output of the frequency discriminator. As a demon- 
stration, we select a typical EEG signal as the input 
waveform to our digital system. We show both the 
original EEG signal and the synthesized EEG signal 
from these two components in Figure 2, Note the high 
fidelity of the reproduction. The filter in the feed- 
back path was chosen to minimize the first order 
difference |[X(N) - X(N-1)||. The EEG data were 
furnished by Dr. J. E. Anliker of NASA-Ames Research 
Center. 
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Figure 1 Block diagram of the digital system 
for resolving the signal into in-phase 
and quadrature components* 



Figure 2 Representation of EEG Signal 

1. x^(N) , the in-phase component* 

2. x 2 (N), the quadrature component* 
3* The original EEG signal. 

4. The reproduced EEG signal 
synthesized from x and x . 

X £ 

5. The reference signal. 
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IN t KM IKK ENT TIME AVKHAtiKS OF KKC DATA 


Jnaes Aniiker, Ph-.D.; David Lai, P h.D.*; 
and Herbert Finger, M.S.E.E. 

NASA-Ames Research Center , Moffett Field 


Tamara Rimmer, B.S 
California 94035 


Coherent time averaging la widely need in tha 
wtudy of «vok»d response » a* « technique for ex- 
trartiog coherent nignalo frura ! nob*- rent noiae. 

11»Jm paper U concerned with tbv problem of o*Ct- 
aj.it Ing, I n «Wl time, tin* signal and noise comp n- 
nerita in coherent time averages of sample functions 
of EEC alpb* activity In the presence, and In the 
absence of coherent trains of photic sclouli. 

The rationale for using a coherent time aver- 
aging algorithm in the study of evoked potentials 
is weli known and seemingly straight-forward . 
events in each sample function that are coherent 
with the time base of the averager C±.e.* bear a 
- consistent temocrai relationship to the trigger) 
v.Lii be enhanced by averaging (or summing) where-^ 
hS incoherent or random events will tend to cancel 
•.Ait. The implication is that when a small signal 
with a lot of noise, the aignal can be 
, x .[ ,1-.^; os* ,:v,*r aging a efficiently large number 
i-« f .uu ;>!« ; on*, t ions. However, in the enfle of elec- 

logical potentials, neither the aignal 
I .,r,*cto: ist'-cs nor the s ignal- to-noise ratio is 
ILwo; consequently, it is not clear how to select 
p^st .'.usetr of sample functions to include in. 
t h« ensemble average. It is also unfortunate that 
.... pri.'o.vb Li^ty that 3 physiological procesa will 
v.'.:.ain stati-iP.ary throughout the averaging time 
cercu^-s as the averaging tint increases. Thus, 
att«r=st tc cuHpeosate for a low signal-to- 
.VV^. 4 -„- £s) iiy incrtuairtg number of sample 

in the av->noi<-. the experimenter may be 

. ; |WM t4 ,,j j,y * , s »n loiitt nmi-Ki.at toiuirlry. Simllur- 
’ friil thvor.'t ft:** l point nl view, lb« 

v-'or imtn Ct; in is double bind: a compromise 
’V^t be "reached between the desire to avoid deci- 
sions -in Insufficient data and the conflict- 

ing desire maintain the shortest possible feed- 
Dack loop. We have sought to minimize these pro 
blems by developing a real time signal- Co-noise 

estimation algorithm. ^ 

We have used a LlNC-S computer for implementa- 
tlon of our estimation scheme and for analyzing 
the coherent changes in the EEC alpha rhythm in 
response to coherent trains of 10 microsec photic 
flashes (Crass PS-2 phocost imulator at intensity 
tit ,) delivered through closed eyelids. The ln *j er “ 
f loH |i period was «el equal to the mean period of 
the autonnswiiH «lj>h.. rhythm as menuured by *uto- 
coi-rvlm hi.u Individual right and 3el t occipital 
sculp oUsirtrodcu were referred to the yoked e«r- 
lobe- electrodes. The cortical potentials were amp 
Ilf led filtered (bandpass- 5 Hz with 24 db per 
octave* roll-off), and converted into digital form 
at * ; msec sampling rate. Each sample function 
was J?*i Tn.sec duration anu was always synchronous 
..j £ -iAsh 'it T . The average amplitude of the 

peaks was measured. This value was treat- 
...Wtho .mixlmiim signal contribution to the aver- 
co hi* expected It the signal were fully co- 
The. dltti repuncy between the coherent pre- 
di-rion m:U the obtained average wn« l rvatmi as 
nt> Ue. The ratio of signal to the combined aig 
nul and noise estimates was re-computed ae each 
n Lw sample function was entered into the average. 


An analog voltage proportional to this ratio was 
Immediately output and displayed nn a storage 
oaci 1 loacope . Figure 1 shows the superimposed 



1 . 0 - 


0 . 5 - 


0 . 0 - 


SUCCESSIVE SAMPLES OFF 


FIGURE 1. Signal /(Signal + Noise) Ration 


successive ratio values for seven 58-semple func- 
tion averages of stimulus-ON (top) and the some 
number of sample functions for the stimulus -OFF 
condition (bottom). The successive ratios for^a 
single cumulative average constitute a sort Oi 
"stabilization pathway*’ for the average. It is 
readily apparent in these records that there is 
greater redundancy as the number of sample func- 
tions in the average increases, afid that there is 
a largo difference between the stimulus-ON and 
ut Imulua-OFF com) 1 ciorin . Consequently, we believe 
that this estimation uch^e hold* coAitlrtoroMr po- 
tential for quantitative, fully automatic ansCs*-. 
Die nta of signal and noise components In this type 
or evoked response. 


The authors wish to acknowledge the technical 
assistance of Mis* Janice McMll lin(NASA/ARC) - 
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A NONLINEAR MODEL OF EEG ENTRAINMENT 
BY PERIODIC PHOTIC STIMULATION 

J. R. Nickolls 

Stanford University 3 Stanford 3 California 
D. c. Lai 

Stanford University and University of Vermont 

J- E. Anliker 

NASA/Ames Research. Center 3 Moffett Field, , California 


The phenomenon that the EEG alpha rhythm synchronizes with periodic 
visual stimuli has been observed since the beginnings of electroen- 
cephalography The brain's entrainment ability depends on several 

factors, including its internal state or condition. Since a quanti- 
tative estimate of .a particular brain state such as alertness has 
many practical applications, a mathematical model of an associated 
input-output relationship is of great use. The quantitative variations of 
the model parameters can be correlated with various brain states using 
empirical data and model simulation. 

We are using a nonlinear oscillator to model a stimulus-response rela- 
tionship of the brain. In particular, we are investigating the entrain- 
ment of the human alpha rhythm by periodic photic stimuli. The model 
elucidates the relationships between the frequency, phase, and amplitude 
of the alpha waveform and the frequency and amplitude of the stimulus. 

This model is clearly applicable to other phenomena such as the flicker 
effect and photically induced epileptic seizures. 

The nonlinear oscillator used in our study is a van der Pol oscillator. 
This has also been suggested by Dewan 2 . The analysis treats en- 
trainment by pulse trains since our experimental stimuli are strobo- •. 
scopic flashes. The theoretical analysis is made with an extension of 
Blaquiere's time-domain method 3 . In this paper, we will report on 
harmonic entrainment, and combined frequency oscillations. Both model- 
simulated results and extensive experimental results from human subjects 
will be described and compared. 

The simulations of the model have been made on a digital computer, and 
the results will be presented for direct comparison with results derived 
from EEG data. The empijfical data were obtained from alert subjects 
photically driven over a wide range of frequencies, and exhibit harmonic. 
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s ubnarmoni c , and superharmonic entrainment, as well as <-nmh-ir,oa -f 
quency oscillations. The model accounts for these critical oh? 

EE^d t tS V E U ft ty haS bee " establishad through extensive empirical"* 
gresf! 0 " "° del and its plications arj i“pro- 
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Abstract 


PREDICTION OF EEG ALPHA WAVEFORMS,, 
USING AN AUTOREGRESSIVE MODEL 

A. Shah, D„ C« Lai, and J. E. Anliker 


An autoregressive process X^ of order p is defined by the 

following difference equation 

X ~ aX 4- a X 4- ... 4* a X 4- e 
t 1 t-1 2 t-2 p t-p t 


where is additive noise and the a r s are the parameters of the 


process. These parameters can be estimated from the time-series by a 
simple non-iterative least -squares method. For the case where the noise 
is Gaussian, these least-squares estimates approximate very closely 
the m ax i mum- likelihood estimates. Such a model is well known for the 
analysis of stationary time-series. In particular, it has been used 
quite extensively for the estimation of EEG power spectra. . However, 
the alpha rhythm of the EEG has slowly changing characteristics which 
can be treated as piece-wise stationary but long samples cannot be 
treated in this manner. We have developed an algorithm which computes 
the least-squares estimates of the parameters of a 7th order process 
from a short sample of the alpha rhythm. The order has been estimated 
from the autocorrelation and partial autocorrelation functions of the 
signal. The algorithm then uses the model for least-squares prediction 
up to some maximum lead time. Probability limits for the forecasts 
(assuming Gaussian noise) are also computed. As more data becomes 
available, the sample window is. moved forward and the parameter estimates 
updated to reflect the changing characteristics of the alpha rhythm. 

The performance of this predictive model is evaluated by the mean-square 
error in the forecasts. The maximum forecast lead time can be varied. We 
have found the model quite useful for prediction as much half an alpha 
cycle (approximately 50 ms.) in advance. 



